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Abstract: Using the Hubeny-Rangamani-Takayanagi (HRT) conjectured formula for 
^ entanglement entropy in the context of the AdS/CFT correspondence with time-dependent 

backgrounds, we investigate the relation between the bulk null energy condition (NEC) 
of the stress-energy tensor with the strong sub-additivity (SSA) property of entangle- 
ment entropy in the boundary theory. In a background that interpolates between an 
AdS to an AdS-Reissner-Nordstrom-type geometry, we find that generically there al- 
ways exists a critical surface beyond which the violation of NEC would naively occur. 
However, the extremal area surfaces that determine the entanglement entropy for the 
boundary theory, can penetrate into this forbidden region only for certain choices for the 
mass and the charge functions in the background. This penetration is then perceived 
as the violation of SSA in the boundary theory. We also find that this happens only 



when the critical surface lies above the apparent horizon, but not otherwise. We conjec- 
ture that SSA, which is thus non-trivially related to NEC, also characterizes the entire 
time-evolution process along which the dual field theory may thermalize. 
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1 Introduction 

Black holes have become, to a modern day's theoretical physicist, an instructive toy 
to play with: the "harmonic oscillator" a la mode. The very nature of black hole 
entropy, which states that the number of degrees of freedom in a theory of quantum 
gravity scales as the area, gave birth to the idea of holography [1, 2]. Subsequently, a 
concrete realization of the holographic principle was conceived under the name of the 
AdS/CFT correspondence [3]. This correspondence, a hitherto unproven 1 but compelling 
conjecture, has also emerged to be a remarkable tool in addressing aspects of strongly 
coupled large N gauge theories. 



1 It is often debated what a proof might mean in this context. In fact, it has been suggested that this 
correspondence should perhaps be viewed as a definition of quantum gravity in AdS-space and not as 
a proposition that is amenable to proof. 
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Within this context, it is possible to realize ideas that are natural in quantum 
information theory: one such example is entanglement entropy. Entanglement entropy, 
defined as the von-Neumann entropy with respect to a reduced density matrix, measures 
the quantum entanglement of a system, and thus becomes an interesting quantity to 
analyze specially for systems described by strongly coupled quantum field theories. For 
large N gauge theories, whose gravity duals are described by Einstein gravity (with a 
negative cosmological constant) in the presence of a suitable matter field, entanglement 
entropy can be computed using the Ryu-Takayanagi (RT) conjectured formula[4] for 
static backgrounds, and later generalized in the Hubeny-Rangamani-Takayanagi (HRT) 
formula[5] for time-dependent backgrounds. The conjectured RT formula has passed 
various non-trivial checks [6, 7] known in quantum information theory and has also found 
numerous intriguing applications; see e.g. [8] for a recent review. 

One important property satisfied by the entanglement entropy is known as the strong 
sub-additivity, henceforth abbreviated as SSA. A quantum system can be described by 
the density matrix, usually denoted by p, which is a self-adjoint, positive semi-definite, 
trace class operator. The entropy of the corresponding system can be described by the 
von-Neumann formula: S = — trfplogp]. 

Let us imagine a quantum field theory on a Lorentzian manifold and further imag- 
ine a Cauchy surface that divides the entire system in two sub-systems: A and A c 
respectively 2 We can now define a "reduced" density matrix for the sub-system A 
by tracing over A c : pA = tr^c [p], and subsequently define a von-Neumann entropy: 
Sa = — tr [pA log pa] as the entanglement entropy. We can now imagine partitioning 
the Hilbert space by more than one Cauchy surfaces. Specifically, if we have three 
sub-systems A±, A 2 and A3, then SSA is defined as 



^A 1 uA 2 + &A 2 uA 3 ~ <->A 2 — &A 1 uA 2 l)A s > , (1-lJ 

Sa 1 ua 2 + Sa 2 ua 3 — Sa 1 — S Aa > . (1.2) 



This property was originally proved in [10, 11], for a recent expository account see 
e.g. [12]. This inequality, that stands as a cornerstone of quantum information theory, 
can be viewed as a crucial ingredient in characterizing the von-Neumann entropy [13, 14]. 
In AdS/CFT correspondence, in a (d + l)-dimensional bulk theory the RT formula 
to compute entanglement entropy of a region A is given by 

S a = ^+n min [ Area (7a)] , (1.3) 



*G N 



(d+i) 



where G N is the bulk Newton's constant, 7^4 denotes the (d — l)-dimensional minimal 
area surface that satisfies d'-fA = dA. For backgrounds with time dependence this 



2 Note that there can be multiple Cauchy surfaces resulting in the same partitioning of the Hilbert 
space. Thus the Hilbert subspace is specified by the (future) Cauchy horizon rather than the Cauchy 
surface itself [9]. 
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proposal is generalized to, via the HRT formula, considering extremal surfaces rather 
than a minimal one. 3 In [6], a simple geometric proof was constructed showing that 
the RT formula obeys the SSA condition, further substantiating the validity of the RT 
formula itself. 

On the other hand, time-dependent backgrounds do provide a more non-trivial check 
of the SSA condition. The prototypical example is the so called AdS-Vaidya background, 
which describes the collapse of a null dust and the formation of a black hole in an asymp- 
totically AdS-space. In the dual field theory this corresponds to a global quench process 
corresponding to the time evolution from a "low temperature" state to a thermalized 
state at a higher temperature. In such a time-dependent background, it was shown in 
[15, 16] that violation of SSA is strongly tied to the violation of null energy condition 
(NEC) from the bulk point of view. 4 For this example, the background is characterized 
by a time-dependent mass function and the NEC imposes a condition on this function. 

Before proceeding further, let us mull over a curious aside. Within AdS/CFT, the 
importance of the NEC has been realized elsewhere: in constructing a monotonically de- 
creasing central charge function along an RG-flow[17-19] for a CFT living in arbitrary 
dimensions. For a (1 + l)-dim CFT, it can be shown that the SSA condition indeed 
implies the Zamolodchikov c-theorem[20]; for more recent developments in higher di- 
mensions see e.g. [21]. Thus, fundamental "inequalities" in a large N gauge theory, 
e.g. a monotonically decreasing central charge along an RG-flow or the SSA condition, 
seem to be stemming from the NEC condition in the bulk description. 

In this article, we intend to sharpen the connection of the SSA condition with the 
NEC condition by studying the formation of a charged black hole in AdS-space. In 
the dual field theory, this will correspond to a global thermalization process in the 
presence of a chemical potential [22, 23]. The corresponding background is a Reissner- 
Nordstrom-Vaidya background in AdS-space, henceforth abbreviated as AdS-RN-Vaidya 
background. This background is characterized by time-dependent mass and charge func- 
tions and the corresponding null energy condition has subtle implications. For a given 
mass and a given charge function, the null energy condition yields a critical surface, 
denoted by z c , that separates the entire spacetime in two regimes: for z < z c , the NEC 
is satisfied and for z > z c it is violated. 5 Therefore, a violation of the NEC depends on 
whether the regime z > z c is accessible to an asymptotic observer. 

It was argued in [25, 26] that for arbitrary 6 mass and charge functions, time-like and 
null geodesies cannot penetrate the critical surface and hence from a gravitational per- 
spective the NEC is protected by having a no-go constraint on these geodesies. However, 



3 In case there are more than one extremal surfaces, one chooses the surface with the minimum area. 

4 A violation of the NEC violates the bound in (1.1), whereas (1.2) remains satisfied. Thus, from a 
holographic perspective, there is a clear distinction between the inequalities in (1.1) and (1.2). 

5 We are working in a coordinate where the boundary of the spacetime is located at z — > 0. 

6 The functions are not completely arbitrary; namely, we need to still impose the same condition on 
the mass function that the NEC imposes in the uncharged AdS-Vaidya case. 
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in AdS/CFT correspondence, space-like geodesies 7 are also relevant since they carry the 
information about non-local operators in the dual field theory such as a 2-point function, 
Wilson loop or entanglement entropy. In this article, we will discuss various examples 
where space-like geodesies can or cannot penetrate this critical surface depending on the 
choices for the mass and the charge functions. In the dual field theory, this penetration 
is perceived as a violation of the SSA condition. Thus, we cannot conspire to have a 
large N gauge theory with certain choices for the mass and the charge functions: entan- 
glement entropy knows it all. However, we will merely discuss generic and instructive 
examples rather than attempting for a general characterization of these functions. 

This article is divided in the following parts: We begin with a short review about 
the SSA and the NEC condition in the AdS-Vaidya background in section 2. Then, we 
discuss in details the AdS-RN-Vaidya background in section 3. We also discuss generic 
examples relating the physics of the NEC condition with the SSA condition based on 
our numerical explorations. We provide examples for asymptotically AdS4 and AdSs- 
backgrounds. Finally, we conclude in section 4. 

2 A brief review 

We begin by briefly reviewing the results that are already known in the literature, 
specially in [15, 16]. 

2.1 Strong subadditivity, concavity and monotone-increasing 

Let us begin by demonstrating the relation of concavity and monotone-increasing with 
the SSA conditions. We will follow closely the discussion in [16]. Let us consider three 
adjacent single intervals Ai, A 2 and A3, whose lengths are denoted by 01, a 2 and 03. By 
symmetry of the construction, S{AA = S(a>i). Now, let us assume that S is a concave 
function. By definition 

S (y Xl + (1 - y)x 2 ) > yS{x x ) + (1 - y)S(x 2 ) , < y < 1 . (2.1) 

Now, let us choose y = 03/(01 + a 3 ) 

xi = a 2 , x 2 = ^a»,, ^ S(a 1 + a 2 )>ys{a 2 ) + {l-y)S l^di) ,(2.2) 

^a; =4> S(a 2 + a 3 ) > (1 -y)s(a 2 ) + yS I ^a; j . (2.3) 

Adding (2.2) and (2.3) we get (1.1). 

On the other hand, the condition of monotone-increasing yields 

5(01 + a 2 ) > S( ai ) , S(a 2 + a 3 ) > S(a 3 ) , (2.4) 



x 2 = a 2 , xi 



7 Through this work, when we mention geodesies, we actually mean extremal area surfaces. 
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adding which we get (1.2). Thus the SSA conditions are equivalent to concavity and 
monotone-increasing conditions on entanglement entropy. 

2.1.1 AdS-Vaidya background and entanglement entropy 

The AdS-Vaidya metric in a (d + 1) spacetime is given by 

L 2 
ds 2 = — [-f(z, v)dv 2 - Idzdv + dx 2 ] , / = 1 - m(v)z d , (2.5) 

which describes the formation of a black holes as a shell of null dust collapses. Here, 
m(v) is a function that interpolates between empty AdS and an AdS-Schwarzschild back- 
ground as a function of v. Also, L is the radius of curvature, x is a (d — l)-dimensional 
vector. We have expressed the above background in the Eddington-Finkelstein coordi- 
nates, where the coordinate v is defined as 

dz 

dv = dt — —, r, (2.6) 

and t denotes the boundary time. In this coordinate system, the boundary is located at 

The energy-momentum tensor that sources this metric has only one non-vanishing 
component: 

T vv = ^z^dMv) . (2.7) 

The null energy condition imposes the following constraint on m(v): 

T^rfn" > =► d v m(v) > , (2.8) 

where n^ is a null vector. 

We now use the HRT formula (1.3) to compute the entanglement entropy for a spatial 
region A. Let's assume A to be a (d — l)-dimensional "rectangle" in the boundary such 
that x 1 G (—£/2,£/2) and x 2 ,...,x d ~ 1 G (0,£j_) at some fixed boundary time t b . The 
HRT prescription establishes that Sa is proportional to the area of the extremal surface 
7a, parametrized by v(x) and z(x), and whose boundary coincides with the boundary 
of A at z = 0. 

Thus, we extremize the area 



/•</2 i . 

Area(7 A ) = L d ~ x V / dx-^y/l - [1 - m(v)z d ]{v') 2 - 2v'z>, ' = d/dx, (2.9) 
J-e/2 z 



n£/2 

-£/2 

where V = £^ 2 ■ We also impose the boundary conditions 

v(-£/2) = v(£/2) = t b , z{-£/2) = z{£/2) = . (2.10) 
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The two equations of motion that follow from (2.9) are 

= [1 - m(v)z d ]v" + z" - dvm{v \ d {v'f - dm{v)z d - l z'v' , (2.11) 

0= zv »-i^l m ( v ) z <i( v >f + (d-l)[( v >) 2 + dv'z'-l] , (2.12) 

and the momentum conservation corresponding to the cyclic coordinate x results in the 
equation 

l-[l-m(z)z d }(v') 2 -2v'z' = (^y {d ~ l \ z , = z(0). (2.13) 

It can be shown that only two of the above three equations are independent. Thus, 
we can solve one of the equations of motion together with the conservation equation 
by imposing appropriate boundary conditions. It is more practical to use the infra-red 
boundary conditions z(0) = z*, v(0) = t>* 8 to solve the equations and, then, read off 
the values of tf, and / from (2.10). Once we solve the system, the extremal area can be 
computed by simplifyng (2.9) using (2.13): 

S(£) = Area( 7 ^) = 2L d ~ l V J <fc^=I) " ( 2A4 ) 

This area contains the usual divergent pieces and we will focus on the finite part only. 

2.1.2 Strong subadditivity and the null energy condition 

In order to illustrate the results found in [15, 16] about the relationship between SSA 
and NEC, we consider two explicit forms of the function m(v) in equation (2.5): 

mi(u) = y(l + tanhf£)) , (2.15) 

m 2 ^) = y(l-tanhf£)) , (2.16) 

which depicted in Fig (1). It is, then, easy to realize that the null energy condition in 
equation (2.8) is obeyed by (2.15) and violated by (2.16). We now specialize in the case 
d = 3 and solve equations (2.12) and (2.13) for both functions m{y) using M = 1.0, 
t>o = 0.01. Then, we plot the entropy function in equation (2.14) for different values of 
the boundary time. It is found that for the case obeying NEC, S(£) is a monotonically 
increasing function that is also concave; whereas for the NEC violating function, S(£) is 
still increasing monotonically but it is not a concave function. These results are shown 
in Figure (2). The fact that S(£) is not a concave function leads to conclude that SSA 
is violated, which establishes a direct connection between SSA and NEC. 



8 Note that these boundary conditions are guaranteed because of the symmetry of our construction 
under x 1 — ► —x 1 . The smoothness of the surface at z* also imposes: z'(0) = and v'(0) = 0. Thus, we 
have sufficient number of boundary conditions altogether. 
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Figure 1. Examples of two functions m(v) where (a) NEC is obeyed, and (b) NEC is violated. 
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Figure 2. Entropy function for the cases where (a) NEC is obeyed, and (b) NEC is violated. 
The different colors correspond to boundary times tb = 0.3 (red), tb = 1.0 (blue), % = 1-5 
(purple), and % = 2.0 (black). Notice that the curves in (b) are not concave and, thus, SSA is 
violated. 



3 The AdS-RN-Vaidya background 

We will now delve into discussing how this connection of SSA and NEC manifests itself 
when there is a non-zero background charge. Here we will flesh out all details in their 
full glory. 



3.1 The bulk action and the backgrounds 

Our initial goal is to write down a metric which describes the formation of a charged 
Reissner-Nordstrom (RN) black hole in AdS-space. One such time-dependent back- 
ground that smoothly interpolates between pure AdS to AdS-RN background is given 
by the so called AdS-RN-Vaidya background. In (d + l)-bulk dimensions, this back- 
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ground is given by 

L 2 
ds 2 = — (-f(z, v)dv 2 - 2dzdv + dx 2 ) , A v = q(v)z d ~ 2 , (3.1) 

fM = l- m{ v)z< + <j0>0l^-\ A=-«^>, (3.2) 

where L is the radius of curvature, z is the AdS-radial coordinate, x is a (d — 1)- 
dimensional vector, A v is the gauge field and m(v) and q(v ) are the mass and the charge 
functions that depend on time. As before, we are working with the Eddington-Finkelstein 
coordinates. The mass and the charge functions, denoted by m(v) and q(v ), are hitherto 
unconstrained. 

The background in (3.1) can be obtained as a solution to the Einstein-Hilbert- 
Maxwell action with a negative cosmological constant coupled to an external source 

S = So + k S ex t , (3-3) 

1 



So 



SnG N 



(d+i) 



l - j d d+l x^~g {R-2K)-- A j d d+1 x^-gF^ 



(3.4) 



where So denotes the Einstein-Hilbert-Maxwell term, S ext denotes the external source 
and k is some coupling. 

The equations of motion resulting from this action are given by 

R,u ~\{R- 2A) 9,u ~ 9 ap F p ,F au + ^g^F^F * = 2 (87rG^ +1) K:) T£ , (3.5) 

d P {V^gg^g^F^) = (svrG^ 1 ^) j; xt , (3.6) 

where the external energy-momentum tensor T®f and the external current J^ xt are 
contained within the action S ex t. More precisely, in order for the equations of motion to 
be satisfied, we have 



2kT- = — ^— - -^z 2d -* q (v)^ S,J VV , (3.7) 



-xt _ / d-l d _ x dm _ d-2 M _ 3 dq 
■>» ~ V 2 dv L 2 q[ >dv 

K J^ t = (d-2)L d - 3 ^ . (3.8) 

In the special case for d = 2, we get 

L 2 
ds 2 = — (—f(z, v)dv 2 — 2dvdz + dx 2 ) , A v = q{v) log z , (3.9) 

z 1 

f(z,v) = l-m(v)z 2 + q -^z 2 \ogz , A = --L . (3.10) 

9 We are considering the case d > 2. The case of d = 2 is somewhat special, which we will briefly 
comment on later. 



The above background is sourced by the following energy-momentum tensor and vector 
current 



2 (dm 2 dq \ 1 dq 

2 (, *T " T^ v W ogz ) ^ ' KJL = LcTv l 



2kT " = o \T7 ~ T-2^)^\ogz) 5, V 5 VV , kJ^ = T ^ . (3.11) 



There is a subtlety in the identification of the source and the VEV in this case and the 
chemical potential should be identified with the sub-leading term rather than the leading 
term of the gauge field as one approaches the boundary. For a detailed discussion on 
this, see [27]. However, we will not discuss this case. 

The energy-momentum tensor presented in equations (3.7) and (3.11) corresponds 
to the energy-momentum tensor of a charged null dust. The easiest way to understand 
this is to note that 

T^-^K, with k 2 = 0, (3.12) 

where we have chosen the vector k^ = 8^ v , which is a lightlike vector. 

Now, NEC — which a reasonable energy-momentum tensor should obey — is given 
by the following inequality: Tj^n^n" > 0, where n^ is lightlike, i.e. n^n^ = 0. There 
are two solutions to the null normal equation n^n^ = 0. Without any loss of generality 
we can write them as 

^ 1)= (0,1,0) , nfa=(l,-±f,8\ , (3.13) 

where denotes the components along the x-directions, which we have chosen to set to 
zero. The null vector n^-, imposes a trivial constraint, and the null vector n^ imposes 

d - 1 rf , dm d — 2 2d _n , .dq ^ n , . tn . .. 

z d 1 - — -z M d q(v)^r > , for d>2, 3.14) 

2 dv L l dv 

^-TU(^log^>0, for d = 2. (3.15) 

dv L 2 dv 

Clearly, NEC is obeyed for all z > z c , where z c denotes the radial position beyond which 
the null energy condition is violated. This critical surface is given by 

j 9 d — 1 L 2 ml , 

* = d=2Ttf> for rf>2 ' (3 ' 16) 

]og z e = ^-^j, for d = 2. (3.17) 

2 qq 

Here ' = d/dv. 

A few comments are in order. Note that for the neutral case, we have q(v) = = 
q'{v) identically. In that case, the critical surface does not exist. In turn, the null energy 
condition then imposes a condition on the mass function: m'(v) > 0. In [15, 16], this 



-9- 



condition was related to the strong sub-additivity property of entanglement entropy, 
which we reviewed in the previous section. In the absence of charge, however, the null 
energy condition seems to be correlated with other simple observations as well. We will 
discuss these momentarily. 

Before proceeding further, let us introduce the apparent horizon for the backgrounds 
in (3.1) and (3.9) following the notations in [28]. The apparent horizon is given by the 
null hypersurface which has vanishing expansion of the outgoing null geodesies. For 
these backgrounds, the tangent vectors to the ingoing and the outgoing null geodesies 
are 

2 2 

£- = -8 z , £ + = -^d v + ^fd z (3.18) 

such that we satisfy 

£_.£_ = o = £+■£+ , £-■£+ = -!. (3.19) 

The codimension 2 spacelike hypersurface, which is orthogonal to the above null geodesies, 
has an area: X = (L/z) 1 . The expansion parameters associated with this hypersurface 
are 

9± = C± log E = t ± d» (log E) , (3.20) 

where C± denotes the Lie derivatives along the null directions £±. The location of the 
apparent horizon is then obtained by solving 6 = 0, where G = 9_9 + . In this particular 
case, the equation 9 = implies f(z,v) = 0. 

In the absence of any charge, we write down a general treatment including the d = 2 
case. The equation for determining the apparent horizon then yields 

1 - m(v)z d = => z^ = m(v)- 1/d . (3.21) 

Here z^ denotes the apparent horizon. Note that in the future infinity, i.e. v — >• oo, the 
apparent horizon coincides with the actual event-horizon. 

Note that, during the time-evolution, a global event-horizon exists in the back- 
ground. This is generated by null geodesies in the background and is the boundary of 
a causal set. Since the background in (3.1) has (d — 1) Killing vectors (d/dx) a , the 
location of the event-horizon is given by a curve z(v ). The null geodesic equation in the 
background (3.1) is given by 

dz eh (v) 1 

d = -7,f( z ch{v) } v) , (3.22) 

where z e h denotes the location of the event-horizon. In the limit v — > +oo, we have 
-2ah = -£eh; however, this is not true in the v — > — oo limit, i.e. the event-horizon lies 
above the apparent horizon. 
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It was argued in [28] that, during a time-evolution, it is the apparent horizon rather 
than the event-horizon that can define a "thermodynamics". Based on an analogy, we 
can define a "temperature function" and an "entropy function" in terms of the apparent 
horizon 



r <") = -££'<*■«> 



= ^mM 1 '" , (3.23) 

•^ah 



S(v) = (V Rd - 1 )m(vp d - 1 ^ d . (3.24) 

Here V^-i denotes the volume of the x-directions. The temperature function is obtained 
by computing the surface gravity at the apparent horizon and the entropy function is 
obtained as the area of the apparent horizon. Clearly, T(v) and S(v) have well-defined 
thermodynamic meaning in the limit v — > +oo. 

Now, taking derivative of these functions with respect to v, we get 

dT{v) dSjv) dm 

— ~ — - — ~ — — > (J , {6.2b) 

dv dv dv 

where we have used the constraint coming from the null energy condition in (2.8). From 
the perspective of the boundary theory, if it makes sense to talk about a "temperature 
function" or an "entropy function" as defined above, the null energy condition implies 
that these must be monotonically increasing. We already remarked that the null energy 
condition was demonstrated in [15, 16] to imply the strong sub-additivity property of 
entanglement entropy. Thus, either all the above observations are physically equivalent 
or we are unable to separate them for the example we are considering here. 

Before proceeding in to the actual computations, let us make some more observations 
here. If we introduce a charge in the system, the null energy condition no longer imposes 
a simple constraint on the mass or the charge function. Instead, it seems to give rise to 
a critical surface z c , above which the energy condition is violated. It was argued in [26] 
that, for such charged backgrounds, null geodesies never penetrate the critical surface 
z c and thus the apparent pathology is not relevant. Within the context of AdS/CFT 
correspondence, spacelike geodesies are also relevant since they contain informations 
about non-local operators, such as 2-point function or the entanglement entropy itself. 
Our goal here will be to analyze further what choices of mass and charge functions 
actually violate the null energy condition and how this is perceived from the perspective 
of the boundary theory as a violation of SSA. For now, we will discuss the case when 
m! > 0, which will smoothly connect to the known results when the charge vanishes. 

3.2 Tests of strong subadditivity 

We now proceed to study entanglement entropy and the SSA inequality in holographic 
theories dual to (d + l)-dimensional AdS-RN-Vaidya spacetimes. As mentioned before, 
d — 2 is somehow special so in this section we will restrict our attention to the cases 
with d > 3. 
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Our starting point is the metric given in (3.1)-(3.2). For simplicity, we consider 
that the region A in the boundary theory is an infinite "rectangular strip" with X\ G 
(—£/2,£/2) and Xj G (—00,00), Vi ^ 1. We will call x± = x given that this is the 
only relevant direction and we will denote the transverse directions collectively as x±. 
According to the covariant prescription for entanglement entropy, we have to find the 
surface 'Ja living in a constant-^ slice that extremizes the proper area functional Area^^). 
This surface is invariant under translation in x±. Thus, without loss of generality, we 
can parameterize it with functions z(x), v(x) and boundary conditions 

z(±£/2)=0 and v(±£/2)=t. (3.26) 

These boundary conditions impose that the boundary of 7^4 coincides with the boundary 
of A along the boundary temporal evolution. The area of this surface is given by the 
following functional 

r i/2 dx 1/9 

S{£) = Area( 7 A) = VL d ~ l / -^ (l - fv' 2 - 2v'z') 1/2 , (3.27) 

J-e/2 z 

where V is the volume that result from integrating over the x± directions. Since there 
is no explicit x-dependence in the Lagrangian, the corresponding conservation equation 
is given by 

/ 7 \ 2(d-l) 

1 - fv' 2 - 2v'z' = (-) , (3.28) 

where z* is defined through z(0) = z*. The two equations of motion obtained by ex- 
tremizing the area functional are 

z +vf + zv — + -v —=0, (3.29) 

az 2 ov 

zv" + (d-l) (v' 2 f + 2z'v' - 1) - -zv' 2 ^f = . (3.30) 

2 az 

In particular, note that the first equation is independent of the dimensions. By taking 
the derivative of the conservation equation (3.28) with respect to x and using one of these 
two equations of motion, one obtains the other one. Thus, it is sufficient to consider 
only (3.28) and e.g. (3.29) to find z(x) and v(x). We, then, solve numerically these two 
equations subject to the boundary conditions 

z(0) = z*, z'(0) = 0, v(0) = v,, t/(0) = 0. (3.31) 

In practice, however, we start the integration at some arbitrarily small x = e to avoid 
possible numerical issues. Also, due to the symmetry of the problem, it is sufficient to 
integrate only for positive values of x. 

So far, z* and v* are two free parameters that generate the numerical solutions 
for z(x) and v(x). The boundary data {£, £{,} can be obtained from these numerical 
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solutions through z(£/2) = z and v(£/2) = t b , where z is a UV cutoff. This cutoff is 
needed because, the area functional (3.27) is divergent and one needs to regularize. The 
divergence comes from the fact that the volume of any asymptotically AdS background 
is infinite and the spatial surface A we are considering reaches the boundary. 

The divergence term can be isolated by studying the same problem in AdS^+i in the 
standard way [29]. Parameterizing with functions x(z) and v(z), it is clear that near the 
boundary x'(z — )■ 0) = 0, v'(z — > 0) = and therefore 



Subtracting this divergence, we obtain the finite term of the area which is the main 
quantity we are interested in, 

S„M = 2VL" £ J(l- fv° - W) 1/2 - — ^ ■ (3-33) 

In order to find numerical solutions to the system (3.28)-(3.29), we employ a "shoot- 
ing" method. First, we give initial values z(0) = z*, v(0) = t>* and integrate until the 
the functions hit the boundary. Once we have the profiles z(x) and v(x), we read 
the boundary values and extract £ and £&. The numerical implementation of (3.33) is 
straightforward. 

In the remaining part of this section, we will consider specific functions for m(v) 
and q(v) in d = 3, 4. As advertised in the introduction, we will find that, although there 
is always a critical surface above which the null energy condition is violated, for some 
appropriate choices of mass and charge functions the extremal surfaces attached to the 
boundary never cross into that region and SSA is satisfied. 

3.3 Thin shells and junction conditions 

Before proceeding to specific examples, let us gain some insight into the properties of 
the critical surface (3.16) let us consider the particular case of a thin null shell located 
at v = which is the boundary of two static spaces. The conditions to join space-like or 
time-like hypersurfaces demand that the two spacetimes induce the same metric on the 
hypersurface and relate the surface stress energy tensor Sab to the jump of the normal 
extrinsic curvature K a b across the hypersurface. The issue is more subtle for null shells 
since the extrinsic curvature no longer carries any transverse geometrical information. 
For null hypersurfaces, the extrinsic curvature is given by tangential derivatives of the 
metric and is, thus, necesarily continuos across the shell and cannot be related to the 
stress energy tensor of the shell S a b- A general formalism applicable to null hypersurfaces 
was developed in [24]. 

For simplicity, consider d = 3. The spacetime metric is given by 

L 2 
ds 2 = -j {~f(z, v)dv 2 - 2dzdv + dx 2 ) . (3.34) 
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Consider two static backgrounds M~ and M + , with mass and charge parameters M$, Qi 
and Mf,Qf respectively. Let M~ and M + be bounded by hypersurfaces X_and X + . 
We glue the two spaces by identifying XL = X + = S. In the present case, we take X 
to be the hypersurface v = 0. We are interested in the case when Mr and A4 + are 
vaccuum solutions with, 

M(v) = M i + Q(v)(M f -M, i ) 1 Q{v) = Qi + Q{v){Q f -Qi)i (3.35) 

f(z,v) = l-M(v)z 3 + ^Q(v) 2 z 4 (3.36) 

= i - (M, + e(«)(M, - m^z 3 + 1 -{q 2 + e( w )(gj - Q 2 ))z A , (3.37) 

where Q(v) is the Heaviside step function. Note that (3.35) allows for the initial and 
final backgrounds to be AdS, AdS-Schwarzchild or AdS-RN. The stress energy tensor 
(3.7) is, 

2 (dM dQ\ 

2K,T fiV = z I — zQ(v)— J b^ v d uv 

= 8{v)z 2 ({M f - M^ - z(Q f - Qi)(Qi + Q(v)(Q f - QO))^^ ■ ( 3 - 38 ) 

Note that (3.38) identically vanishes in M~ and M + and is non-zero only at v — 0. 
This discontinuity comes with a sound physical interpretation; it is associated with the 
presence of a thin distribution of matter at v = 0. The only non-zero component of the 
surface stress energy tensor is T zz , 

2kT£ u = 5(v)z 2 ((M f - Mi) - z(Q f - Q t )(Qi + Q(v)(Q f - Q l )))k^k v (3.39) 

= 5(v)a(z)5» z 5 uz . (3.40) 

Since there is no rest frame for a null shell, we cannot formally identify a as the surface 
density. However, a can be used to determine the results of measurements by any spec- 
ified observer. This involves introducing an arbitrary congruence of timelike geodesies 
intersecting X associated with the different families of observers making measurements 
on the shell. An observer with four- velocity u a = dx a /d,T, will measure an energy den- 
sity associated with the shell T^ u u^u v = 5{v)a{z){k^u^) 2 . Thus, the arbitrariness of the 
choice of congruence is limited to an overall factor and the quantity cr(z) is independent 
of this choice. It is in this sense that we interpret a(z) as the shell's surface density. 

From (3.16), we see that for this type of backgrounds the critical surface is located 
at v = and 

7 M f~ M > _ 2(M f -M t ) 

c (Q f -Qi)(Qi + e(o)(Qf-Qi)) Qj-Ql ' K ' ' 

where we have used O(0) = 1/2 as is conventional in distribution theory. Evaluating a 
at the critical surface we obtain, 

ff(Zr) = . 
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Thus, the critical surface is the locus where the shell's surface density becomes zero. 

Our main interest is to study extremal spacelike surfaces in the backgrounds de- 
scribed above. From the point of view of the spacelike surface, there is a discontinuity in 
the dz/dv when the surface crosses the shell. This can easily be seen from the equations 
of motion. As in the previous section, consider a rectangular strip in the boundary the- 
ory. Extremizing the area functional (3.27) we obtain the equations of motion (3.29-3.30) 
which for d = 3 read, 

zv " + Av i z > + (2/ - ||>'2 = , (3.42) 

z" + fv" + fz'v' + lfv ,2 = 0. (3.43) 

dz 2 dv 

The z(x) and v(x) coordinates are continuos across the shell while / and -£ present a 
finite jump. From (3.42) we see that z', v' and v" remain finite. On the other hand, 

£ = z 3 [(M f - M t ) - z{Q f - Qi){Qi + B(v)(Q f - Q t ))}5(v) 

= -f(z,v)S(v) (3.44) 

diverges at v — 0. Thus, the behavior of (3.43) across the shell is, 

z" ~ -hf( ZjV )S(v)v' 2 (3.45) 



and z' has a finite jump. Indeed, 



[z'V - (z') 



o+ ,0+ z n 1 



o- Jo 



z"dx= —dv f(z,0)v' (3.46) 



v' 2 



and the jump of ^ across the shell is, 



*®-@) + - ©"->■«> 



where + and — refer to quantities evaluated in Ai + and Ai respectively. Using (3.41) 

dz 
dv 



we find that the critical surface corresponds to the surface at which the jump on 0: 



vanishes, 



A (l 



~f(z c ,0) = 0. (3.47) 



Intuitively this makes sense; we know that at z c the mass of the shell goes to zero, the 

dz 
dv ' 



shell has disappeared and there is no reason for a jump in d: 



-15- 



3.4 Examples in d = 3 

Let us begin our discussion for the (3 + 1) -dimensional bulk theory, where the dual 
field theory is a (2 + 1) -dimensional conformal field theory in the presence of a chemical 
potential. Presumably the corresponding UV-completion is given by an if? -reduction 
of 11-dimensional supergravity, which leads to an 5*0(8) gauged supergravity in (3 + 
l)-dimensions. Therefore, the boundary theory should correspond to an ABJM-like, 
i.e. Chern-Simons-matter theory in the presence of a chemical potential. 10 

3.4.1 Backgrounds that respect SSA 

We will begin with the examples that do not violate SSA. A good starting point is to 
use the functions analyzed in [23] that were used to address scaling properties of the 
thermalization time with respect to the temperature and the chemical potential of the 
thermalized state. These are: 

m(v) = - ( 1 + tanh ( 

v ; 2 V V0.01 

It is clear from Figure (3) that there is no change in concavity of the entanglement 



and 



q(v) =0.9m(V) 2/3 



(3.48) 
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Figure 3. Left panel: m(v) (black) and q(v) (red). Right panel: S(£), for % 
0.5 (blue), 1 (purple), 1.5 (magenta), 2 (green), 3 (orange) and 5 (red). 



0.01 (black), 



entropy function as £ grows for a given boundary time £&. 

According to equation (3.16), there exists a critical surface in the bulk geometry 
beyond which NEC is violated. It is thus instructive to analyze whether the space- 
like geodesies, which eventually determine the entanglement entropy, can penetrate this 
critical surface or not. In Figure (4), we show a representative family (characterized by 
the length of the entangling region) of geodesies corresponding to function (3.48). We 
also display the critical surface (3.16) and the location of the apparent horizon. 



10 We should note that while this is a very plausible scenario, we are making an assumption that the 
Vaidya-type backgrounds can be embedded within gauged supergravity consistently at least in some 
well-defined limit. 
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Figure 4. Profiles of a family of geodesies when SSA is obeyed. The dashed gray line represents 
the apparent horizon at v = and the red line represents the critical surface. This family is 
parametrized by £, the length of the entangling region at the boundary. Note that right panel 
shows that the geodesies do not intersect the critical surface as explained in text. 



Before going further we should note that the critical surface, in the thin-shell limit, 
exists within a tiny 11 region around v — 0. Also, we note that the coordinate v evolves 
along each geodesies independently and at v = these geodesies cross the shell, which in 
the thin-shell limit coincides with the apparent horizon. Figure (4) therefore compares 
the location of the critical surface to the location of the geodesies at v = 0. 

There are two key features that stand out from Figure (4): First, the critical surface 
lies above the apparent horizon 12 ; in other words, it is cloaked by the apparent horizon. 
Second, although the space-like surfaces cross the apparent horizon, they do not probe 
the forbidden region beyond the critical surface irrespective of how large £ becomes. At 
this point we would like to stress that the above observations seem very robust against 
a substantial amount of test cases. One might imagine designing a situation where the 
critical surface comes arbitrarily close to the apparent horizon, thus encouraging the 
geodesies to cross it. However, we have verified that this does not seem to happen. 

Let us now illustrate a couple of more examples where the SSA condition is satisfied. 

and q(v) = — ( 1 - tanh ( 777T7 ) ) . (3.49) 



m[v) 



1 + - ( 1 + tanh , 

2 V V0.01 



m(v ) = - ( 1 + tanh ( 

K J 2 V V0.01 



and 



q(v) 



0.9 



( tanh ( 

V V0.01 



0.01 



tanh 



1 



0.01 



(3.50) 

The choice functions are given in equations (3.49) and (3.50) and the corresponding 
figures are shown in Figure (5) and (6) respectively. 

Before concluding this section, let us offer some remarks. In view of (3.25), it can be 
verified that the condition dT(v)/dv > is not satisfied for all v with the choices made in 



12 



'The width of this region is characterized by the shell thickness parameter vq. 

1. 6. Z c S> Z a f l . 
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Figure 5. Left panel: m{y) (black) and q{v) (red) given in (3.49). Right panel: S(£), for t\, 
0.01 (black), 0.5 (blue), 1 (purple), 1.5 (magenta), 2 (green), 3 (orange) and 5 (red). 
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Figure 6. Left panel: m(v ) (black) and q(v) (red) given in (3.50). Right panel: S(£), for tf, = 
0.01 (black), 0.5 (blue), 1 (purple), 1.5 (magenta), 2 (green), 3 (orange) and 5 (red). 

(3.50). Thus the SSA condition is an independent constraint which is not related to the 
rate of change of surface gravity at the apparent horizon in the dynamical geometry. As 
far as the corresponding physical processes are concerned, the choices in (3.48) represents 
a situation in the dual field theory, where both temperature and chemical potential are 
increasing from a "low value" to a higher non-zero value. On physical grounds, this is 
perhaps the most "reasonable" process. 

The choices in (3.49) takes a low temperature, high chemical potential initial state 
to a high temperature, low chemical potential final state. Finally, the choices in (3.50) 
are rather exotic, which takes a low temperature, vanishing chemical potential initial 
state to a high temperature vanishing chemical potential final state; but does not obey 
dT(v)/dv > 0, Vf. As far as NEC or SSA is considered, there is nothing preventing 
these two choices; however, whether they are realizable as solutions of gravity with a 
reasonable matter field is an issue we will not address here. 
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3.4.2 Backgrounds that violate SSA 

Now let us illustrate a few examples where SSA is violated. One such choice is: 
0.05 ^ . / u \\ . 0.9 



m\v) 



0.95 



1 + tanh 



0.01 



and 



q(v) 



1 + tanh 



0.01 



(3.51) 

The corresponding illustration is shown in Figure (7). The corresponding geodesic is 
shown in Figure (8). Once again, the critical surface exists only around a small neigh- 
bourhood of v = and the geodesies reach the shell at v = 0. Note here that both the 
features alluded to in the previous subsection are gone: First, the critical surface lies 
outside the apparent horizon at v = 0; second, the minimal area surface brings news 
from the forbidden region in the bulk by probing the region beyond the critical surface. 
This is perceived as the violation of SSA condition in the boundary theory. At this 
point we emphasize that these observations seem rather generic and hence we will not 
pictorially illustrate a similar behaviour of the geodesies for other representative cases, 
whenever SSA is violated. 
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Figure 7. Left panel: m(v) (black) and q(v) (red) given in (3.51). Right panel: S(£), for % 
0.01 (black), 0.5 (blue), 1 (purple), 1.5 (magenta), 2 (green), 3 (orange) and 5 (red). 



Let us now take a second example where the violation of SSA is observed: 

v-Y 



m(v) = - ( 1 + tanh ( 

v ; 2 V V0.01 



and 



q(v) 



0.9 
~2 



tanh 



0.01 



(3.52) 



In this case also the minimal area surface probes the region beyond the critical surface. 
The corresponding plots showing the violation of SSA are presented in Figure (9). 

Once again, a few comments are in order. First, it can be checked that for the 
choices in (3.52) we will get dT(v)/dv > and dS(v)/dv > and still a violation of 
SSA. Thus, the SSA condition is indeed independent of these. Intuitively, the process 
in (3.52) is not very meaningful since it seems to allow for the black hole to accumulate 
charges keeping it's mass fixed. On the other hand, there is no such a priori objection 
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Figure 8. Profile of geodesic when SSA is violated corresponding to the choice (3.51). The 
dashed line represents the apparent horizon or the position of the shell at v = and the red 
line depicts the critical surface at v = 0. Clearly the geodesic probes the region behind the 
critical surface at v = 0. 
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Figure 9. Left panel: m(v) (black) and q(v) (red) given in (3.52). Right panel: S(£), for % = 
0.01 (black), 0.5 (blue), 1 (purple), 1.5 (magenta), 2 (green), 3 (orange) and 5 (red). 

to the process in (3.51), and still SSA is violated. This indicates that SSA is a rather 
non-trivial condition on the allowed trajectories a thermalization process might take for 
a given field theory. 

3.5 Examples in d = 4 

We will now move up by one dimension and consider a (3 + l)-dimensional conformal 
field theory. The dual geometry will correspond to the formation of a charged black 
hole in five-dimensional AdS-space. In this case, the possible UV-completion will be 
given by an S^-reduction (or a reduction on a Sasaki-Einstein five-manifold) of type IIB 
supergravity truncated to the M = 2 sector with an U(l) 3 symmetry (or at least one 
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£7(1)). Thus the dual field theory is presumably a cousin of the prototype N = 4 super 
Yang-Mills (SYM) theory in the presence of a chemical potential. 13 

3.5.1 Backgrounds that respect SSA 

We will discuss similar processes as in the case for d = 3. The qualitative features 
are very similar here, and hence we will limit ourselves in terms of the details. The 
analogous choices that preserve SSA are pictorially represented in Figure (10) and (11). 
Note that, Figure (10) corresponds to the choices used in [23] to analyze the scaling of 
the thermalization time with temperature and chemical potential of the system. 
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Figure 10. Left panel: m(v) (black) and q(v) (red). Right panel: S(£), for t^ = 0.01 (black), 
0.5 (blue), 1 (purple), 1.5 (magenta), 2 (green), 3 (orange) and 5 (red). 
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Figure 11. The mass and charge functions for which SSA is satisfied: m(v) (black) and q(v) 
(red). 

13 Once again we note that this claim needs to be rigorously demonstrated, which we will not attempt 
here. 
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3.5.2 Backgrounds that violate SSA 

Once more, we will keep our discussion very brief and present the examples that violate 
SSA. These are shown in Figure (12). The choices are analogous to the ones made in 
(3.51) and (3.52). This physics can also be observed by studying the geodesies, which 
penetrate the critical surface whenever there is a violation of SSA but not otherwise. 
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Figure 12. The mass and charge functions for which SSA violation occurs: m{v) (black) and 
q(v) (red). 



4 Discussion and conclusions 

In this article we have explored and demonstrated an interesting pictorial realization of 
the strong sub-additivity condition in terms of the bulk gravitational description. In 
the presence of charge, the dynamical evolution from a low temperature, low chemical 
potential pure state to a thermalized state with a non-zero value of the chemical potential 
does indeed sharpen the connection between the bulk null energy condition constraint 
and the strong sub-additivity of entanglement entropy in the boundary field theory, 
which was alluded to in [15, 16]. 

Our investigations suggest that the dual field theory disallows specific choices of the 
mass and the charge functions for which it is possible to penetrate the critical surface. 
However, as we have learned now, the nature of the violation of the NEC depends on 
the class of examples we choose; such as the details which have qualitatively distinct be- 
haviours for the charged case and the uncharged one. Generally, the NEC is an algebraic 
constraint on the bulk energy-momentum tensor, which — should a general result exist 
- may correspond to an algebraic constraint in the boundary theory as well. Strong 
sub-additivity depends crucially on the concavity property of the entropy function and 
thus it is an intriguing possibility to consider establishing a direct equivalence between 
the bulk null energy condition and the concavity property of the entropy function at 
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the boundary. See [30] for some recent progress towards a proof, however, it does not 
necessarily apply for backgrounds where a black hole eventually forms. 

Coming back to our case, it is perhaps surprising how SSA can be obeyed for some 
examples, specially since the critical surface always exists for any generic mass and 
charge functions. Unlike the time-independent cases, where no extremal surface can 
penetrate the black hole event horizon, the Vaidya backgrounds give rise to an apparent 
horizon that can be penetrated by a space-like surface. The critical surface may lie 
above or below this apparent horizon. If it is below, the violation of NEC is imminent. 
However, there is no a priori criterion that prohibits the minimal surface to penetrate 
the critical surface when it is cloaked by the apparent horizon. Nevertheless, this is 
what we observe. It will thus be interesting to unveil whether this holds true for just 
the specific examples we considered, or perhaps a general feature. 

Nonetheless, based on our limited observations we can venture a naive character- 
ization for the choice functions. Let us assume that m'(v) > 0. This is required by 
continuity with the results known in q(v) —¥ limit. Given this, we can characterize 
a charge function q(v) to be good if z c > z^ for all times and bad if z c < z ah for any 
time. Of course, we also need to impose a constraint on the maximum magnitude of the 
mass and the charge functions in order to avoid the naked. Such a characterization, at 
present, is only a plausibility. 

It is intriguing that the SSA condition seems to constrain the global time-evolution 
process but does not say anything about the initial conditions. In general, it is possible 
that given an arbitrary but reasonable initial condition the dual field theory undergoes 
time-evolution, but never thermalizes or obtains a steady-state phase. Such a process, 
once obtained by solving Einstein gravity with a reasonable matter field in the bulk, will 
surely preserve NEC and hence SSA. It will thus be interesting to investigate whether 
the SSA condition plays a similar role for more conventional systems rather then large 
N gauge theories having a gravity dual. We hope to address these issues in details in 
future. 
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